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ABSTRACT 

We present an analysis of NICMOS photometry and low-resolution grism 
spectroscopy of low-mass stars and sub-stellar objects in the young star-forming 
region NGC 1333. Our goal is to constrain the ratio of low-mass stars to sub- 
stellar objects down to 20 Mj up in the cluster as well as constrain the cluster 
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IMF down to 30 Mj up in combination with a previous survey of NGC 1333 by 
Wilking et al. Our survey covers 4 fields of 51.2" x 51.2", centered on brown 
dwarf candidates previously identified in Wilking et al. We extend previous 
work based on the use of a water vapor index for spectral typing to wavelengths 
accessible with NICMOS on the HST. Spectral types were derived for the 14 
brightest objects in our fields, ranging from <M0 - M8, which at the age of 
the cluster (0.3 Myr) corresponds to a range in mass of >0.25 - 0.02 M . In 
addition to the spectra, we present an analysis of the color-magnitude diagram 
using pre-main sequence evolutionary models of D'Antona & Mazzitelli. Using 
an extinction-limited sample, we derive the ratio of low-mass stars to brown 
dwarfs. Comparisons of the observed ratio to that expected from the field IMF 
of Chabrier indicate that the two results are consistent. We combine our data 
with that of Wilking et al. to compute the ratio of intermediate-mass stars (0.1 

- 1.0 M ) to low-mass objects (0.03 - 0.1 M Q ) in the cluster. We also report the 
discovery of a faint companion to the previously confirmed brown dwarf ASR 
28, as well as a possible outflow surrounding ASR 16. If the faint companion is 
confirmed as a cluster member, it would have a mass of ~ 5 Mj up (mass ratio 
0.15) at a projected distance of 350 AU, similar to 2MASS 1207-3923 B. 

Subject headings: stars: pre-main-sequence, brown dwarfs - stars: mass function 

- stars: formation - infrared: stars - ISM: individual (NGC 1333) 



1. Introduction 



The shape of the Initial Mass Function (IMF) and its connection to the initial physical 
conditions in molecular clouds remains one of the fundamental questions in star formation. 
More specifically: Does the young cluster IMF mimic the integrated field star IMF even to 
very low masses? Is there a low-mass cut-off to the sub-stellar mass function? 

The IMF over the full range of st e llar m asses has been extensivel y studied in t h e past 
deca des starting wit h iMiller fc Scald (Il979l ). as well as updates by iKroupa et al.l (119931 ) 
and IChabrierl (120031 ). It is generally accepted that young clusters exhibit a Salpeter-like 
IMF above 1 M , flattening out towards low-mass stars. Studies of the IMF have been 
concentrated on young clusters, because they offer several beneficial characteristics. Their 
populations are less likely to have undergone significant dynamical mass segregation com- 
pared to older clusters, which means that a small field can yield a sample representative 
of the whole cluster. In addition, young low mass stars are still located significantly above 
the main sequence and are thus several times brighter than their main sequence counter- 
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parts. Lastly, young clusters are fairly compact, meaning that they occupy a small area in 
the sky as well as often being situated in front of a large amount of extinction. This mini- 
mizes both foreground contamination by field stars and background contamination by giants. 



Recent studies of galactic young clusters have been extended down t o the hydrogen 



burnin g limit to a distance of 1 kpc (ILuhmanl (119981 ): iHillenbrandl (119971): ICarpenter et al. 
jl997h ). These studies suggest that the IMF is universal above 0.1 M e jMever et al.lboooh . 



The shape of the IMF in the sub-stellar r egime is less well constra i ned a n d we are only 



now b eginning to probe it in great detail ([Hillenbrand fc Carpenter! (120001 ); iBriceno et al. 



(120021 )). Several clusters have been studied down to 30 M j up in both photome tric and 
spectroscopic s tudie s, among them th e Orion Nebula Cluster (ISlesnick et al.l (120041)), IC348 



( Naiita et al. 



200 



1 



Luhman et al.l (j2003al )). the Taurus star for ming reg ion 



Luhman 



k^Uranan er, al ] t ha Chanson T star forming re gl on Cta3653)) and 

the Pleiades JMoraux et all (booah ). 



Searching for IMF variations at the lowest masses can tell us about possible low-mass 
cut-offs in the IMF, which m ight be asso ciated with the minimu m mass for opac i ty lim ited- 
frag mentation 10 Mjup) (IReesl 119761 ) . Evidence presented by IBriceno et al. rl2002h sec 
also: ILuhmanl (j2004bl )) showed the Taurus star-forming region to exhibit a dearth of brown 
dwarfs, by as much as a factor of 2, when compared to the denser region of the Orion Nebula 
Cluster. However, newer studies, covering a larger area of Tauru s suggest that this deficiency 
is not as marked as previously reported (ILuhman et al.l 120061 ). Nonetheless the low-mass 
IMF remains a strong candidate for variations with cluster characteristics. This suggests a 
possible connection between stellar density and the shape of the low-mass IMF. NGC 1333 
is a young star-forming region intermediate in density between Taurus and Orion, making it 
an ideal test bed to search for variations in the low-mass IMF and explore differences based 
on star-forming environment. 



The NGC 1333 reflection nebula is part o 



estimated distance of 300 pc (Ide Zeeuw et al.l (119991 ); iBelikov et al.l (120021 ) ). The proximity 



the P erseus molecu l ar clo ud complex at an 



of the cluster allows us to study the IMF dow n to very low ma s ses. I t was first identified as 
a star- forming region in Herbig Sz Rao ( 19721 ) (see also: Herbig ( 1974 ): Strom et al. ( 1976T)). 



Infrar ed surveys have sin ce revealed a large pop ulation of low-mass objects ( lAspin et al. 
(119941 ): lLada et al.l (119961 ): IWilking et al.l (120041 ) ). These studies have characterized the 



young stellar population over large areas of the cloud, but have not fully probed the very 
lowest masses because of a lack of sensitivity. 

We have performed a new, deep (J < 21 m ) near-infrared survey over a limited area in 
NGC 1333 using NICMOS/HST. This study is complete to fainter magnitudes than previous 
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studies, thus allowing the characterization of a significant extinction-limited sample of brown 
dwarfs. This enables us to study the population of very low-mass objects in detail. We 
obtained photometry as well as infrared spectroscopy. The spectra enable estimates of stellar 
masses and ages for individual objects by comparing their positions in the H-R Diagram with 
PMS evolutionary models. The spectra further help us adopt models to interpret our flux- 
limited survey by constraining a mass-luminosity relation appropriate for the cluster. We 
construct an extinction limited sample to explore the shape of the mass function in the 
cluster. 

The paper is structured as follows. Section 2 describes the observations and data reduc- 
tion, followed by a description of the photometric and spectroscopic analysis in Section 3. 
Section 4 details the results, and places them in context with recent literature. A summary 
and conclusions are presented in Section 5. 



2. Observations and Data Reduction 

We obtained images of the young cluster NGC 1333 using camera 3 of NICMOS (NIC3) 
as part of the HST program 9846. Six fields, each 51.2" x 51.2" with a plate scale of 
0.2"pixel _1 , were obtained between 2004 Jan 15 and 2005 Aug 11 using 12 HST orbits. In 
addition to the 6 NGC 1333 fields we also obtained observations for 7 additional objects. 
Six of these are previously confirmed brown dwarfs from older star-forming regions and the 
field. We obtained these observations in order to explore the surface gravity dependence of 
our spectral typing technique. The coordinates for each of the fields in addition to the total 
exposure times are listed in Table 1. For the position of our fields with respect to molecular 



line maps of the region please see Figure 1 in IWilking et al.l (120041 ) . Each field was observed 
in F110W, F160W and the grism G141, which cover 0.8 - 1.4 jjm (broadband), 1.4 - 1.8 
/im (broadband) and 1.1 - 1.9 /jm (dispersed), respectively. F110W is roughly equivalent 
to J-band and F160W is roughly equivalent to H-band. The fields were observed in a 2 
x 2 dither pattern, with an offset of 1.75 pixels. The observations were subdivided into 
visits, with each visit covering one field in G141, F160W and F110W in that order. Since 
the grism G141 is slitless, stars that are close in either the x or the y direction can have 
overlapping spectra for each individual visit. To minimize this overlap we observed each 
photometric and spectroscopic field in 3 visits, at roll angles offset 30 - 75° and 105 - 150° 
between subsequent visits. We were unable to obtain 3 separate roll angles for the fields N2 
and N7 causing their final photometry to be less deep than that of the other 4 fields. For this 
reason they are only used in the spectroscopy section and excluded from the photometric 
analysis. All 4 of these fields lie in the southern half of the cluster and the complete area 
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covered by our photomet r ic sur vey i s 2.9 arcmin 2 . Com pared to the combined areas of the 
surveys by IWilking et al.l (120041 ) and lAspin et al.l (Il994j ) , which surveyed both the north and 
south components of the cluster and covered 79 arcmin 2 and 81 arcmin 2 respectively, our 
photometric survey thus covers ~ 2% of the area of NGC 1333. Total integration times for 
the fields with 3 roll angles, including dithers and multiple roll angles were 383 s for F160W, 
766 s for F110W and 4608 s for G141 per field. 



2.1. Photometry 



Data reduction was carried out using a combination o f IRAF and custom IDL routines. 
The methods described below closely follow those used in |W. Liu et al.l (120031 ). The images 
were dark and sky subtracted using combined dark plus sky frames cre ated with the rou- 
tine NICSKYDARK in the NICRED package for IRAF jMcLeodlE997l ). Cosmic rays and 
bad pixels were located and removed using the routine FULLFITBAM, by searching for 
discontinuities in the flux in each pixel over time. After dark and sky subtraction and cos- 
mic ray reduction, there did not appear to be any bias offsets or "pedestal effects" between 
the different q uadrants of the NICM OS chip, which is a common problem experienced with 
NICMOS (e.g. IW. Liu et all {^J). This is because our fields were not very crowded and 



did not show significant nebulosity, as the pedestal is also dependent on the total charge 
of the quadrant. The only field that shows any nebulosity is S2. Finally, the images were 
flatfielded using the routine NICFLATTEN and the appropriate epoch on-orbit flatfields 
from the Space Telescope Science Cente r Institute (STScI ) website. The dithers were then 
combined using the IDL software IDP3 (ILytle et al.l 119991 ) . aligning the different frames by 
their World Coordinate System as well as resampling the images by a factor of 2 using a 
bi-cubic sine to increase the resolution of the images. The three roll angles were combined 
using IDP3, first rotating the images to align in angle and then shifting the images to correct 
for offsets in x and y. For both of these routines, flux was conserved. 

All sources presented in the Color-Magnitude Diagram (CMD) (Sec. 3.3) were detected 
in both the F160W and F110W filters using the IRAF routine DAOFIND with a 10 c 
detection threshold. Even with the low amount of nebulosity every image still had some 
spurious detections which were removed by visual inspection. We took advantage of the 
three roll angles for each field to help identify false detections. If a source appeared in the 
same position for each roll angle, it was assumed to be genuine. False detections were located 
at their position on the chip for only one roll angle. 

Since the fields were not crowded, the photometry was performed using the IRAF rou- 
tine APPHOT. The optimal aperture size was calculated to be 8 pixels in radius with the 
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sky background measured using an annulus between 10 and 12 pixels around each object. 
Aperture corrections were calculated out to 25 pixels using several bright, non-saturated 
stars and applied to the photometry. The aperture corrections derived were 0.054 ± 0.008 
mag for mllO and 0.073 ± 0.008 mag for ml60, where mllO and ml60 refer to the mag- 
nitudes associated with F110W and F160W, respectively. Errors in the photometry were 
computed using PHOT and agreed with errors computed by comparing photometry derived 
from each separate roll angle. Two objects, ASR 9 and ASR 28, showed close companions. 
For these an aperture radius of 3 pixels was used with appropriate aperture corrections of 
0.32 ± 0.03 mag for mllO and 0.38 ± 0.02 mag for ml60 and a sky annulus of 20 to 24 pixels 
so as to exclude flux from the companion. In addition, one object, ASR 16, appears to have 
extended nebulosity associated with it. Because of this, photometry was extracted for this 
object in the same manner as the binaries to minimize contamination from the nebulosity. 
This source is further discussed in section 4.2. A total of 25 unique sources were detected 
and are listed in Table 2. 

NICMOS magnitudes were calibrated using zero-points of 1775 and 1093 Jy and were 
calibrated to the Vega system using the conversions of 2.873 x 10~ 6 and 2.776 x 10 -6 Jy 
ADU -1 s _1 for F110W and F160W respectively. Detection limits were assessed using ar- 
tificial star tests. Artificial stars were added to the images in 0.5 magnitude steps using a 
PSF derived from bright, non-saturated sources. The recovery fraction of these stars was 
then computed using the same detection technique as described above. The 90% complete- 
ness limits were found to be ml60 ~ 20. 5 m and mllO ~ 21. m . The NICMOS magnitudes 
ml 10 and ml60 were transformed to the CIT system, J and H, by comparing magnitudes for 
sources in the survey which also had 2MASS photometry (a total of 12 obj ects). The 2 MAS S 



photometry was converted to the CIT system using the transformations of lCarpenterl (120001 ) . 
A linear regression was then performed to obtain color transformations between the HST 
magnitudes and the CIT system. The following color corrections were determined: 

H = ml60 + (0.189 ± 0.030) + (0.120 ± 0.025)(mll0 - ml60) 
J-H = (0.132 ± 0.073) + (0.760 ± 0.060)(mll0 - ml60) 

These color transformations were used to transform the NICMOS magnitudes to CIT with 
the J magnitudes calculated as J = (J-H) + H. 

Relative astrometry for each object was determined using the World Coordinate System 
in the he ader of the NICMOS images together with the XY2SKY routine in the WCSTools 



package (IMinkl |2002| ). The derived astrometry was compared with previously determined 
coordinates from 2MASS and showed no systematic offset. All coordinates presented in 
Table 3 are estimated to be accurate to < 1.5". 
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2.2. Spectroscopy 

The spectroscopic observations were performed with the grism G141, which is centered 
at 1.4 /iin, with a wavelength coverage between 1.2 /zm - 1.9 /zm and a resolution of ~ 200 per 
pixel. This covers the water-band feature at 1.4 /mi which is expected in the atmospheres 
of cool objects and is difficult to observe from the ground, making HST ideal for these 
observations. To reduce these data, the images were first run through the HST pipelines 
CalnicA and CalnicB. These pipelines perform the dark subtraction, using artificial dark 
frames from the STScI website, as well as cosmic ray reduction. The pipelines were not used 
to co-add either the dithers or the roll angles, because this did not improve the quality of 
the spectra over extracting them from individual images. This means that for each object 
in the NGC 1333 sample there can be a maximum of 12 extracted spectra. The spectra 



were extracted using a custom IDL routine NICMOSlook (jPirzkal et al.lll998l ). which was 
designed to deal specifically with NICMOS grism data. NICMOSlook reduces spectra in a 
standard way, by tracing the spectrum across the chip and summing the flux in the spatial 
dimension as well as subtracting a background region set by the user at each wavelength. 
To extract spectra, NICMOSlook needs a spectroscopic as well as a photometric image, to 
perform the wavelength calibration. The photometric image gives the true location of the 
object's position, which then gives the zero point for the wavelength of the spectrum. The 
photometric image used was the F160W image because its central wavelength is closer to the 
central wavelength of G141. The spectra were extracted with an extraction width of roughly 2 
x FWHM of the spectra, in general ~ 4 pixels. The background size was varied depending on 
the crowding of the field. NICMOSlook also flatfields the spectra during the extraction. This 
is necessary because the quantum efficiency of the detector changes with both wavelength and 
position. For this reason the spectra cannot be flatfielded before extraction. NICMOSlook 
flatfields the spectra using a set of narrowband flatfield images, which are chosen according to 
the date of the observations, from which the software constructs a 3-dimensional calibration 
datacube. 

Signal-to-noise ratios (SNR) were estimated for all extracted spectra in two ways. First, 
a 3rd order polynomial was fitted to the combined average spectrum in the spectral region 
between 1.5 - 1.8 jum and the RMS of the difference between the observed continuum and 
the fit per pixel calculated. Second, a SNR was calculated for each individual pixel in the 
spectrum, by estimating the signal as the mean of each pixel and the noise as the error in 
the mean from the ensemble of individual spectra for each source. Since there are only four 
independent extractions for some of the spectra, the error in the mean is not necessarily a 
high fidelity estimate for the noise in each pixel. This means that the SNR estimated in this 
way may overestimate the real SNR of the spectra. The SNR of the spectrum is then the 
average of the SNR for all the pixels. 
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2.3. Spectral Standards 

To help calibrate our method for spectral typing we obtained spectra of 7 objects with 
previously known spectral types (see Table 1), six of which are brown dwarfs. These objects 
are older than the NGC 1333 sources and were observed to explore the surface gravity 
dependence in our spectral typing technique. Of the 7 objects, only 4 had high enough 
SNR to be included in the spectral sequence. They span a spectral range from M8 to L7. 
In addition we extracted spectra f rom archival HST grism data (Program IDs #7322 and 



#7830; see also lNajita et al.l (120001 )) of old low-mass field dwarfs with known spectral types. 
They span a spectral range from MO - L5. The spectra are plotted in Figures 1 and 2. All 
standards used in the spectroscopic analysis are listed in Table 3. 



2.4. NGC 1333 Spectra 

Spectra were extracted for a total of 14 objects in NGC 1333 from 6 fields. Depending 
on the position of the spectra on the chip, such as proximity to other objects and position 
with respect to the edge of the chip, between 4 and 12 spectra were extracted for each source. 
The spectra were then combined in a sigma-clipped way, excluding pixels which deviated by 
more than three sigma. For the NGC 1333 objects the SNRs ranged from 30 to « 100 and 
the different SNR estimates agreed well with each other. We compared these values with the 
expected SNR from the STScI exposure time calculator. In general these theoretical SNR 
tended to be higher than those estimated from the real spectra. We chose to adopt the SNR 
estimates from the polynomial fits. Spectra of all NGC 1333 objects are plotted in Figure 3. 



3. Analysis 
3.1. Spectroscopic Analysis 



Spectral typing was performed using a water vapor absorption index centered on the 
water band at 1.4 ^m. This band is sensitive to temperature in co ol stars later t h an M O 
(IJones et al.l Il994f) but not strongly dependent on surface gravity (IGorlova et al.l (120031); 
Wilking et alii 2004)). The met hod described here is similar to the one detailed in Wilking et al 
jl999h and lWilking et al.l ( 20041). We defi ne a Q index that is independent of reddening, using 



the reddening law of 

Q = (F1/F2)(F3/F2) - 567 



Cohen et all f]198lh : 
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where Fl, F2 and F3 are the averaged fluxes in the bands 1.30 - 1.35 /im, 1.40 - 1.45 /im and 
1.65 - 1.70 fim respectively. The index is independent of reddening because it includes only 
a ratio of flux bands which have been scaled by an exponent directly related to the applica- 
ble reddening law. This accounts for any inherent slope in the spectrum due to interstellar 
reddening and corrects for it. To optimize our spectral typing technique we chose among 
several different combinations of flux bands and determined which fit of Q versus spectral 
type for our spectral standards showed the lowest combined error in slope and intercept and 
at the same time showed the strongest evolution with spectral type. This index was chosen 
because of its tight fit and its strong evolution with spectral type, as it directly measures 
the evolution in the depth of the water band feature at 1.4 /im. The evolution of Q with 
spectral type for the NGC 1333 objects is shown in Figure 3. The uncertainty in Q was 
estimated by calculating Q for each individual spectrum of an object and then deriving the 
error in the mean of those values from the multiple spectra which were available for each 
object. To calibrate the relation between Q and spectral type, a weighted linear fit of Q 
vs. spectral types for the old field dwarf standards and the young confirmed brown dwarfs 
was performed. Together these standards have spectral types ranging between MO - L7 (see 
Table 3). Spectra with types earlier than MO were excluded because they show no water 
band absorption. This fit gives the relation: 

MV(subclass) = (-5.04 ± 0.31)Q + (-2.22 ± 0.47) 

The fit, together with the Q values for the field dwarfs, is plotted in Fig. 4. 

We tried to assess the surface gravity dependence of the Q index by deriving a fit of Q 
vs. spectral type for the young brown dwarfs with known spectral types observed as part 
of our program. Field stars in general have log(g) = 5.0 - 5.5, while younge r PMS stars 



have lower surface gravities between log(g) = 3.0 - 4.2 (IGorlova et al.l 120031 ). Thus it is 
important to explore whether our spectral typing technique is valid across a range of log(g). 
We performed a fit using only the young confirmed brown dwarfs and compared it to the fit 
derived using the field dwarf standards. Comparing these two fits, we estimate that assigning 
spectral types using the old M dwarfs might underestimate their temperature by ~ 200 K 
(too late by one subcla ss). This means the t emp eratures we der i ve are lower limits, similar 



to previous findings by IGorlova et al.l (120031 ) and IWilking et al.l (120041 ) . We only have high 



SNR spectra for 3 young brown dwarfs however, and so these are preliminary results at best. 

Table 4 shows the derived spectral types for the NGC 1333 objects. The spectral 
types derived for the NGC 1333 objects range from <M0 to M8. Stars with no water band 
absorption features were assigned a spectral type of <M0. Uncertainties in Q were between 
±0.01 - 0.17, which leads to uncertainties in spectral types of approximately one subclass, 
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comparable to the potential error in our calibration mentioned above. From this relation we 
classify seven objects with spectral types of M6-M8 which are likely brown dwarfs. Three of 
the objects for which we have derived spectra, ASR 24, ASR 1 5 and ASR 17, have previous 
spectral classifications using K-band H 2 absorption bands (jWilking et ajj 120041 ) . Both 
derived spectral types are listed in Tab le 4 and agree to wit hin the errors. While we did not 
detect cooler young brown dwarfs than lWilking et al.1 (12004 ) . we were able to observe similar 
objects with greater obscuration. Alternative methods of the spectros c opic r eduction as well 
as a method of spectral typing are also addressed in detail in iFanettil (120051 ). 



3.2. H-R Diagram 



To estimate the age of the cluster we compared the positions of the NGC 1333 PMS 
objects in the H-R Diagram to theoretical evolutionary tracks by ID'Antona fc Mazzitelli 
( 1l997;1998[ ) (hereafter DM98) as appropriate for the median cluster age characterizing the 
low mass population. First, we converted spectral typ es derived f o r the NG C 1333 objects 



to tem peratures using empirical relations derived from lDahn et al.1 (120021 ) and lLeggett et al. 



( 119961 ). The typical uncertainties in temperatures due to uncertainties in spectral types were 
« 200K. Luminosities were derived from the de- reddened J-band magnitudes. To determine 
Mj we first calculated vis ual extinctions using the standard relation A t) =9.09x[(J-H) o6s - 
(J-H)J (jCohen et al. 1981 ) where (J-H) G is the intrinsic color for a given spectral type. The 
values of A„ derived for NGC 1333 sources lie between and 14 mag, with 6 sources having an 
extinction higher than 5 mag. Values for log(L^/L Q ) were calculated from the de-reddened 
J-band magnitudes, after bolometric correction, as l ogfL^/Lfn) = 1 .89- 0.4 M feoi . We adopt 
intrinsic colors and bolometric corrections following IWilking et al.l (120041 ) . Typical errors in 
log(L fcoi /L ) are ± 0.2 dex and are derived based on errors in the J-band magnitude, the 
distance modulus and the bolometric correction. The H-R Diagram is presented in Figure 5. 
The luminosities and temperatures for the NGC 1333 objects with spectral types are listed 
in Table 4. The positions of the brown dwarfs are compared to the theoretical tracks and 
isochrones from the models of DM98. In ad dition to our sample, objects in NGC 1333 for 
which unique spectral types were derived in IWilking et al.1 (120041 ) are plotted on the H-R 
diagram. From these objects we estimated the cluster to be younger than 1 Myr and chose the 
DM98 0.3 Myr tracks as appropriate for the median cluster age. It is noteworthy that these 
ages are relative ages as absolute ages are model dependent. The objects in the HR Diagram 
show a scatter around the median age with an overall trend to younger ages for lower mass 
objects. The scatter is due in part to errors in the luminosities of the objects introduced by 
errors in the photometry and dereddening. The overall trend towards younger ages for lower 
mass objects is an effect often observed in HR Diagrams of young stellar clusters. This is 
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most likely an effect of imperfect PMS models and not a real effect of an age-mass relation 
for the cluster. The median age of 0.3 Myr for our H-R diagram corresponds to the break 
point of the mass bins for our analysis (0.1 M ), thus ensuring that we are using the correct 
mass-luminosity relationship at that mass. 

Four objects exhibit unusual positions in the H-R Diagram, with ages much older than 1 
Myr, appearing under-luminous compared to their PMS counterparts. Two of these objects 
have spectral types earlier than MO, meaning they show no water band absorption feature 
in their spectra. These two objects could be background stars. The other two objects, both 
of spectral type M4 also a ppear fainter than ex pected, but could still be cluster members. 
Notice also 3 objects from IWilking et al.l ( 2004 ) which lie below the expecte d locus for the 
majo rity of the PMS objects. These objects are specifically addressed in IWilking et al. 
(120041 ). but the following reasoning applies to them as well. The sources could have an 



infrared excess, causing a dilution of the water absorption bands, leading us to assign a 
spectral type that is too early. They could also exhibit unresolved scattered light due to a 
surrounding disk, which would result in bluer colors and could cause us to underestimate 
extinction a nd thus lumin o sity. There is evidence that disks around young brown dwarfs are 
ubiquitous. |M. Liu et al.l (120031 ) for example found that 77% of a sample of young brown 
dwarfs in IC348 and Taurus showed evidence of a circumstellar disk. To establish whether 
any of these objects show signs of a disk and thus confirming their youth, we attempt to 
assess whet her they show an inf rared excess in the next section. This method is described 
in detail in lWilking et al.l (119991 ). Since we cannot rule out either solution at this point it is 
impossible to assign a mass to the two M4 sources or an age if they are not cluster members. 
For these reasons all four sources have been excluded from the mass function analysis. It 
is possible that other objects in our sample also possess circumstellar disks. However, for 



most objects we can expect th e effect of the disks on spectral typing to be small (IMeyer et al. 



(119971 ) ; lMuzerolle et al.l (120031 ) ). Even if objects show a strong infrared excess only those with 
special viewing geometry would be dominated by scattered light leading us to underestimate 
the reddening. Further observations will be required to test whether these objects suffer 
from the effects mentioned. 



3.3. K-band excess 

For NGC 1333 objects with 2MASS photometry we have estimated the amount of 
excess emission at K, defined as ^k—^Kex/^K 1 -, where Fxex is the flux contribution at K 



X r fc = F Kex /F K = [(i +rA )(io[( ff -^)-(«- A ') -o.065A„]/2.5)_ 1 ]^ We have assume( j t ] aat r .= r/l =0. 
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due to circum s tellar material and the inherent stellar flux at K (jWilking et al.l (120041 ); 



Wilking et al.l (119991 ) ) . Values of r*. are listed in Table 4. These estimates are lower limits 



since we have assumed no excess at J- or H-band. Two objects, ASR25 and ASR105, show 
moderate K-band excesses. ASR25 has a peculiar position on the H-R Diagram as mentioned 
above. The measured K-band excess makes it a likely cluster member and it may indicate 
that we have assigned it a spectral type that is too early, which would make it a brown 
dwarf, and that it may possess a disk. Most of the other objects show little or no excess 
emission. It is possible that the lowest mass objects possess disks which are too cool to cause 
an excess i n the K-band as inner disk temperature is expected to be correlated with stellar 
luminosity ( Pascucci et al. ( 20031 )). 



3.4. Color-Magnitude Diagram and Extinction-Limited Sample 



In Figure 6 we present a CMD for the cluster. Also shown is the 0.3 Myr isochrone 
of DM98, which we adopted based on the H-R diagram discussed above. The isochrone 
was convert ed to th e CMD by using a set of colors and bolometric correctio ns, with colors 
taken from ILeggett Jl992h. whi l e bol ometric corrections were adopted from ILeggett et al.l 



( 119961 ). similar to lWilking et all (119991 ). All but 3 out of 25 objects in the CMD appear to lie 
below the hydrogen burning limit for a cluster age of 0.3 Myr. This number changes to 5 if 
we instead assume a cluster age of 1 Myr. The photometry, isochrone and reddening vector 
( iCohen et al.lll98ll ) are all presented in the CIT system and the completeness limit discussed 
above is also shown. For the age derived above we define an unbiased sample of low-mass 
objects in order to constrain the IMF for the cluster. We begin by creating an extinction- 
limited subsample, to ensure that we are sampling the stellar population uniformly and not 
over-representing more luminous (massive) deeply embedded objects. The extinction limits 
used are A^ < 21.1 mag and < 18.3 mag for the 0.3 Myr and 1 Myr DM98 isochrones 
respectively, corresponding to a mass range of 0.02 M - 0.1 M . For both isochrones this 
led to a total sample of 13 (17-4) objects, excluding the 4 objects mentioned in section 3.2. 

We are sensitive to a mass range lower than 20 Mj up , albeit for a smaller range of 
extinctions. However, the DM98 tracks do not extend below this mass limit. Thus no 
objects lower than 20 Mj„ „ were considered for the cluster mass function. This means that 
compared to the survey of IWilking et al.l (120041 ) we probe to a higher A„ but not to a lower 
mass range. There are six objects lying in the region below 20 Mj up for an assumed age of 
0.3 Myr, two of them very close to the 90% completeness limit. One of them, ASR 28b, is 
a suspected companion to a brown dwarf primary and will be dealt with in more detail in 
Section 4.3. Another object has been previously classified as ASR 16 and appears to have 



13 



extende d nebulosity a s socia ted wit h it. There are track s available that extend below 20 Mj up , 
such as iBaraffe et al.l (119981 ) and iBurrows et al.l (119971 ). However the NGC 1333 objects lie 
above the youngest tracks available for both sets of models (1 Myr) in the H-R Diagram. 
Thus these tracks are most likely not appropriate for NGC 1333, making the DM98 tracks 
the most appropriate choice for NGC 1333, with a median age of 0.3 Myr. Typical errors in 
assig ning masses from different tracks are approximately a factor of 2 (IWhite fc Hillenbrand 
20041 ) although this is only well tested down to ~ 0.3 M . 



It is also worth noting that while the rapid drop in the number of objects below 20 
Mj up looks striking, a fairly large area in the CMD in that region is covered by a small 
range in mass name ly be tween 5-20 Mj„. p , w ith the masses defined according to tracks by 



Baraffe et al.l (119981 ) and lBurrows et al.l (119971 ). It is however curious that the objects we do 
detect seem clustered towards the lower end of this mass range, which is close to our 90% 
completion limit. We detect no objects between 10 - 20 Mj up with low A„. Due to our small 
sample size and survey area this does not allow us to rule out that such objects exist in the 
cloud. It is however suggestive that NGC 1333 may lack objects below 20 Mj up . 

It is possible that at least some of the objects are background contaminants. We have 
tried to asses s the amount of b ackground contamination by adjusting the number counts 
measured by iLucas et al.l (120031 ) in the Hubble Deep Field South (HDF-S) using F160W. 
Using their estimates for number counts one expects 1.7 +0.4/-1.1 objects between 19 m - 
21 m in our survey. However, the HDF-S points out of the Galaxy and thus this number 
should be a lower limit since NGC 1333 lies closer to the Galactic plane. On the other 
hand our NGC 1333 fields lie in front of a large amount of extinction, which will reduce the 
number of contaminants. 



4. Discussion 
4.1. Cluster IMF 

To constrain the slope of the IMF down to the lowest masses in the cluster we derive the 
ratio of very low-mass stars (0.076 - 0.1 M Q ) to brown dwarfs (0.02 - 0.076 M ) in our survey. 
We calculate this ratio to be R = 3/10 = 0.30 ±0.20 assuming a 0.3 Myr isochrone, with the 
errors computed due to Poisson statistics. For a 1 Myr isochrone this ratio increases to R = 
5 /8 = 0.625 ± 0.3 56. Both cluster ages give ratios consistent with having been drawn from 



a IChabrierl (120031 ) system IMF, which gives a most likely ratio (mode) of R = 0.30. This 
corresponds to a slope of dN/dM oc M _L1 over this range. For the older isochrone there are 
a larger number of low-mass stars but still more brown dwarfs due to the hydrogen burning 
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limit shifting to a fainter magnitude for a higher cluster age. Both ratios are lower limits, 
since the low-mass bin has objects for which we do not have spectra and thus contamination 
by field stars may lead us to overestimate the number of brown dwarfs relative to low mass 
stars in the sample. Above 0.076 M all objects but one have assigned spectral types which 
allow us to assess their cluster membership. Th e object without a spectral type is ASR9a 
which is a previously identified cluster member (lAspin et al.lll994l ) which we have resolved 
to be a binary. Below 0.076 M Q however, apart from ASR 9b, there are four objects without 
spectral types, which could be possible field stars. Thus these ratios imply upper limits for 
the slope of the cluster mass spectrum below 0.1 M of a < 1.1 and a < 0.55 for 0.3 Myr 
and 1 Myr isochrones respectively, where dN/dM oc M~ a and a = 2.35 is the Salpeter slope. 



Where does our result stand with regard to comparable clusters and studies? IWilking et al. 



( 120041 ) have explored the mass spectrum of NGC 1333 using an extinction limited sample 
down to 0.04 M Q finding the ratio of sub-stellar (0.04 - 0.1 M ) to stellar objects (0.1 - 1 
M Q ) to be Rss = 1 11 +0.8/-0.4. They used this to estimate an upper limit for the slope of 
the lower end of the mass function of a < 1.6. This c ompares wel l with slopes in the solar 
neig hborhood be l ow th e hydrogen-burning limit, which iReid et al.l (119991 ) found to be 1 < a 
< 2. lAllen et al.l (120051 ) have further attempted to constrain the shape of the field star IMF 
below 0.1 M and found -0.6 < a < 0.6 with 60% confidence and a best fit of a = 0.3. 



To derive a broader estimate of t 



re clu ster mass function we attempted to combine 



this survey with that of IWilking et al.l (120041 ) to make the cluster mass function more easily 
compara ble to ratios publishe d for other regions. We combined all objects between 0.1-1 
M from lwilking et all J2004f ) with an extinction limit of A v < 12.8 m with the HST objects 
between 0.03 - 0.1 M with an extinction limit of A„ < 21. l m , scaling the number of HST 
objects by the ratios of the survey areas, which is 26.3, as we ll as the extinc t ions, which is 
0.61. We chose not to include the objects below 0.1 M from IWilking et alj (120041 ) as that 
survey was estimated to have a large field star contamination below 0.1 M . We compute 
this ratio as Ri = N(0.1 - 1 M )/N(O.O3 - 0.1 M ) = 0.14 ± 0.047 with errors again due 
to Poisson statistics, using the same cluster age and models for both surveys. This ratio is 
abnormally low when compared to other published results in this mass range even within the 
errors. We have recomputed this ratio for 5 different young star clusters based on previously 
published data sets and results range between « 3 - 6. All surveys were either extinction 
limited surveys or were normalized to make them as close to extinction limited surveys as 
possible. See Table 5 for computed results and errors. 

We then tried to compensate for possible differences in stellar density between the two 
surveys by introducing a normalization factor in an overlapping mass range between 0.07 - 
0.1 M . This normalization factor was computed to be 3.5 ± 2.0. This increases the ratio to 
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R.2 = 0.64 ± 0.38, which is higher but still significantly lower than the result for any other 
cluster. This survey was centered on known brown dwarf candidates. Could this bias our 
result towards an artificially high number of brown dwarfs? Each of our fields was centered 
on one brown dwarf candidate, with S3A and S3B including the same candidate but offset 
from each other. If we remove these three candidates from our ratio, the normalized ratio 
increases to R3 = 0.83 ± 0.41, which is consistent with the ratio previously determined. 
Thus the inclusion of these objects does not fundamentally change our results. 



Due to the expected field star contamination in I Wilking et al.l (I2004J ) both ratios are 
lower limits. However, if the observed density of brown dwarfs is constant over the molecular 
core previously surveyed at near-infrared wavelengths, this would imply that brown dwarfs 
outnumber low-mass stars in the cloud. This may be some indication that brown dwarfs 
form in a clustered environment within the cloud or are spatially segregated from stars due 
to dynamical evolution. 



4.2. Pre-Main Sequence Activity 



Getman et all (120021 ) have conducted a Chandra X-Ray study of young stellar objects in 
NGC 1333. Two objects in our study, ASR8 and ASR24, show signs of x-ray activity. Both 
are spectroscopically confirmed brown dwarfs and have virtually no extinction. This means 
that 2 of our 8 (excluding the 4 sources discussion in Sec. 3.2) or 25% of our spectroscopically 
confirmed brown dwarfs are detected in the X-Ray. However, most of the other spectroscopic 
brown dwarfs in our s urvey have higher ex tinctions, decreasing the chance of detecting 
existing x-ray activity (IPreibisch et al.ll2005l ). 



One of the objects in our survey, ASR 16, appears to have extended nebulosity asso- 
ciated with it, which is detected only in the F160W image. We were not able to obtain a 
spectrum for the object or the associated nebulosity. Extended emission in the infra-red can 
be associated with emission-line Herbig-Haro objects or scattered light envelopes. Since the 
nebulosity is seen only in the H -band and not at J, contrary to expectations from scattered 
light models (IGomez et al.lll997l ) we explore the possibility of an Herbig-Haro outflow associ- 
ated with this brown dwarf candidate in NGC 1333. Since the nebulosity emission is bright in 
the F160W filter which covers 1.4 - 1.8 /im it seems that the emission is likely due to [Fe II] at 
1.644 /im. Such outflow s are thought to be associated with accretion onto the central object. 
Muzerolle et al.l (120051 ) find signs of accretion are common in the spectra of brown d warfs. 
Outflows have also recently been detected around brown dwarfs (IWhelan et al.ll2005l ). We 
tentatively identify this as a candidate outflow around the brown dwarf ASR 16. 
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4.3. Binaries 

There are two sources in our fields which have possible close companions. One was pre- 
viously identified as ASR 9, which is now resolved into ASR9a and ASR9b, with a separation 
of 0.8". The other is ASR28a, which we have determined has a close companion ASR 28b. 
ASR28a is one of our brown dwarf candidates with a spectral type of M7. ASR28a and 
ASR28b have a separation of 1.2" (« 350 AU at the distance of NGC 1333) and ASR28b is 
4.2 magnitudes fainter in H. If ASR28b is indeed a cluster member its colors and apparent 
J-band magnitude (assuming the sam e A„ as ASR 28a) suggest a temperatu re of & 1800 K 



(spec tral type L4 (IDahn et al.l 120021 )) and a corresponding mass of 5 M Jup (IBurrows et al. 
19971 ) (with an approximate mass ratio of 0.15). If confirmed this object would be very sim- 
ilar to the reported planetary mass companion of 2MASS 1207-3923 (hereafter 2MASS1207). 
ASR28b is too faint to extract a spectrum from this survey and J- and H-band photometry 
are not enough to rule out the possibility that it is a background object earlier than M4. To 
ascertain that it is indeed a cluster member we hope to obtain a J- and H-band spectrum of 
ASR28b in the near future. 



5. Summary and Conclusions 



We present NICMOS photometry and spectroscopy f or 6 re gions in NGC 1333 centered 
on previously detected brown dwarfs from IWilking et all (120041) . We also a nalyze spectra of 
brown dwarf standards from this program and that of iNajita et al.l (120001 ). Our results are 
as follows: 



We detect a total of 27 sources at H and 25 sources at J down to Mj < 21 m and M# 
< 20. 5 m over a region of 4 x 51.2" x 51.2" square arcseconds. 



Following IWilking et al.l (120041 ) and IWilking et al.l (119991 ) we develop a reddening- 
independent water band index to estimate spectral types within one subclass from 
SNR > 30 NICMOS G141 spectra. With this method we derive spectral types for 
14 sources in NGC 1333 from <M0 - M8. Seven objects have spectral types > M6 
suggestin g they are young bro wn dwarfs. While we were sensitive to lower mass ob- 
jects than lWilking et al.l (120041 ). we detected only objects of comparable mass but with 
much higher extinction. 

Using the CMD and an assumed age of 0.3 Myr for the cluster we define an extinction 
limited sample of 13 sources between 0.02 - 0.1 M with A„ < 21. l m . Although we 
are sensitive to lower A^ objects between 10 - 20 Mj up , none were detected. Though 
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we are limited by the small number of objects in our survey, this may be indicative of 
a lack of objects below 20 Mj up in the cloud. 

• We compute the ratio of low-mass stars to brown dwarfs in the cluster as R = N(0.076 
- 0.1 M )/(O.O2 - 0.076 M m ) = 0.30 ± 0.20, which is consistent with having been drawn 



from a field-star IMF (jChabrierll2003l ) . We also compute the ratio of sta r s to lo w-mass 
objects in the cluster in conjunction with the survey by IWilking et all (120041 ) as R= 
N(0.1 - 1M )/(O.O3 - 0.1M Q ) = 0.64 ± 0.38, which lies below published results for 
other regions. This may indicate that brown dwarfs form in a segregated environment 
compared to stars in this cluster. 

• Our survey includes several unusual objects, including a possible companion to a spec- 
troscopically confirmed brown dwarf ASR 28, which if it is a cluster member would 
have a mass of 5 Mj up . In addition one of the brown dwarfs (ASR 16) might have an 
outflow associated with it, which could be confirmed through follow-up spectroscopy. 
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Fig. 1.— NICMOS G141 spectra of field dwarf standards from MO - M6. 
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Fig. 2. — NICMOS G141 spectra of field dwarf standards as well as younger brown dwarfs 
confirmed through visible spectroscopy between M6.5 - L7. The younger brown dwarfs are 
denoted by an asterisk. 
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Fig. 3. — NICMOS G141 spectra for candidate brown dwarfs located in NGC 1333, arranged 
from earliest to latest spectral type. 
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Fig. 4. — Plot of Q vs. spectral type showing the fit derived for the old field dwarfs and 
young brown dwarfs. 
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Fig. 5. — H-R Diagram of the NGC 1333 o bjects from ou r surv ey overplotted on DM98 
isochrones and tracks as well as sources from IWilking et al.l (12004] ) . For sources in common 
between the two studies we defer to this study. Four objects show unusual positions in the 
H-R Diagram. They are commented upon in detail in section 3.2. 
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Fig. 6. — NICMOS photometry of s ources detected in NGC 1333 tra nsformed into the 
CIT system together with a 0.3 Myr ( iD'Antona &: Mazzitellilll997; 19981 ) isochrone and the 
reddening vector also in the CIT system. The 90% completeness limit, derived from artificial 
star tests, is also shown. The objects with derived spectral types are marked. 
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Table 1. Log of Observations 



Obs. Field 


R.A. (J2000.0) 


Dec. (J2000.0) 


F110 exp. time 


F160 exp. time 


SI 


03:28:56.41 


+31:15:33.1 


766s 


383s 


S2 


03:29:03.16 


+31:16:58.7 


766s 


383s 


S3A 


03:29:11.97 


+31:16:58.2 


766s 


383s 


S3B 


03:29:12.38 


+31:17:25.2 


766s 


383s 


N2 


03:29:11.65 


+31:23:19.5 


510s 


256s 


N7 


03:29:04.1 


+31:25:29.5 


255s 


128s 


Roque 25 


03:48:30.6 


+22:44:50 


256s 


128s 


Roque 33 


03:48:49.0 


+24:20:25 


256s 


128s 


Teide 1 


13:05:40.18 


-25:41:06.0 


256s 


128s 


J0205.5-1159 


02:05:29.40 


-11:59:29.7 


256s 


128s 


Sigma Ori 60 


05:39:37.62 


-02:30:45.64 


256s 


128s 


Sigma Ori 47/27 


05:38:16.00 


-02:40:23.80 


256s 


128s 



Table 2. Photometry of sources detected in J- or H-band 



Object ID R.A. (J2000.0) Dec. (J2000.0) J Jerr H Herr ASR ID a Notes 



Possible Outflow 



Sl-1 


03:28:54. 


.938 


+31:15:29. 


.65 


16.09 


0.14 


14. 


.96 


0. 


.05 


109 


Sl-2 


03:28:55. 


.334 


+31:15:53. 


,83 


N/A 


N/A 


20. 


.00 





,23 


N/A 


Sl-3 


03:28:56. 


,020 


+31:15:49. 


.28 


20.62 


0.25 


19 


.05 


0. 


,10 


16 


Sl-4 


03:28:56. 


.959 


+31:15:50. 


93 


15.10 


0.13 


14. 


.06 





,04 


15 


Sl-5 


03:28:57 


,162 


+31:15:35. 


.19 


15.28 


0.15 


13 


.95 


0. 


.05 


17 


Sl-6 


03:28:57. 


.388 


+31:15:36. 


.72 


17.83 


0.18 


16 


.30 


0. 


.07 


N/A 


Sl-7 


03:28:57. 


.990 


+31:15:50. 


36 


N/A 


N/A 


21. 


03 





.39 


N/A 


Sl-8 


03:28:58. 


.058 


+31:15:48. 


.65 


19.24 


0.32 


16 


.67 


0. 


.10 


N/A 


S2-1 


03:29:01 


.673 


+31:17:12. 


.25 


20.44 


0.31 


18. 


.96 


0. 


.16 


N/A 


S2-2 


03:29:01 


.876 


+31:16:43. 


.98 


21.79 


0.62 


19 


.79 





.24 


N/A 


S2-3 


03:29:01 


,929 


+31:16:53. 


.22 


16.77 


0.15 


15. 


.47 





,06 


11 


S2-4 


03:29:02 


,654 


+31:17:14. 


.84 


22.44 


0.66 


21 


.45 


0. 


,48 


N/A 


S2-5 


03:29:04 


,077 


+31:16:59. 


.71 


20.95 


0.39 


19 


.36 





,19 


N/A 


S2-6 


03:29:04 


.114 


+31:17:07. 


.46 


13.25 


0.10 


12 


.67 


0. 


.03 


8 


S2-7 


03:29:04 


.179 


+31:16:50. 


.93 


18.74 


0.28 


16 


.56 


0. 


.09 


9b 


S2-8 


03:29:04 


.229 


+31:16:50. 


60 


17.15 


0.24 


14. 


.80 





.09 


9a 


S2-9 


03:29:04 


.713 


+31:16:59. 


03 


15.00 


0.17 


13 


.44 


0. 


.06 


105 


S2-10 


03:29:05 


.413 


+31:17:02. 


.57 


20.56 


0.41 


18. 


.33 





.14 


N/A 


S3B-6 


03:29:10 


.716 


+31:17:21. 


.34 


21.36 


0.44 


20 


.04 





.25 


N/A 


S3A-1 


03:29:10 


.787 


+31:16:43. 


.99 


16.07 


0.17 


14. 


.55 





,06 


23 


S3A-2 


03:29:11 


,272 


+31:17:18. 


,85 


13.95 


0.10 


13 


.30 


0. 


.04 


24 


S3B-3 


03:29:11 


.353 


+31:17:17. 


.52 


13.99 


0.09 


13 


.31 


0. 


.03 


24 


S3A-3 


03:29:12 


,923 


+31:17:08, 


.34 


15.94 


0.11 


15. 


.38 


0. 


.04 


25 



Both in frame S3A & S3B 
Both in frame S3A & S3B 
Both in frame S3A & S3B 



Table 2 — Continued 



Object ID 


R.A. (J2000.0) 


Dec. (J2000.0) 


J 


Jerr 


H 


Herr 


ASR ID a 


Notes 


S3B-4 


03:29:13.008 


+31:17:07.01 


16.02 


0.15 


15.40 


0.04 


25 


Both in frame S3A & S3B 


S3A-4 


03:29:13.014 


+31:16:41.12 


18.19 


0.22 


16.30 


0.08 


26 




S3B-7 


03:29:13.034 


+31:17:37.34 


20.39 


0.32 


18.80 


0.13 


28b 




S3B-1 


03:29:13.078 


+31:17:38.24 


15.24 


0.10 


14.62 


0.04 


28a 




S3A-5 


03:29:13.269 


+31:17:17.35 


20.21 


0.41 


17.54 


0.12 


27 


Both in frame S3A & S3B 


S3A-6 


03:29:13.283 


+31:16:45.91 


23.37 


2.38 


19.47 


0.37 


N/A 




S3B-5 


03:29:13.351 


+31:17:16.01 


20.14 


0.39 


17.55 


0.12 


27 


Both in frame S3A & S3B 


S3B-2 


03:29:13.655 


+31:17:43.48 


16.14 


0.20 


14.19 


0.07 


29 




S3A-7 


03:29:14.285 


+31:17:00.58 


21.05 


0.36 


19.37 


0.19 


N/A 





a ASR refers to lAspin et all (J199J) 
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Table 3. Spectral Standards 



ID 


R.A. (J2000.0) 


Dec. (J2000.0) 


Spectral Type a 


Q b 


Prop ID 


GL 328 


08:55:07.62 


+01:32:47.4 


M0 


0.87 


7322 


GL 908 


23:49:12.53 


+02:24:04.4 


Ml 


0.98 


7830 


GJ 623 


16:24:09.32 


+48:21:10.5 


M2.5 


1.25 ± 0.01 


7322 


GL 388 


10:19:36.50 


+ 19:52:10.6 


M3 


1.38 


7322 


GL 569a 


15:43:02.12 


+26:16:36.0 


M3 


1.50 


7322 


GJ 748 


16:24:09.32 


+48:21:10.5 


M3.5 


1.15 


7830 


GL 896a 


23:31:52.18 


+ 19:56:14.1 


M3.5 


1.46 


7830 


GL 213 


05:42:09.27 


+ 12:29:21.6 


M4 


1.56 


7322 


GJ 699 


17:57:48.50 


+04:41:36.2 


M4 


1.16 ± 0.01 


7322 


GJ 473a 


12:33:16.3 


+09:01:26 


M4 


1.26 


7830 


GL 83.1 


16:36:21.45 


-02:19:28.5 


M4.5 


1.56 


7322 


GJ 1245a 


19:53:54.48 


+44:24:53.3 


M5.5 


1.23 ± 0.01 


7830 


GL 406 


10:56:28.99 


+07:00:52.0 


M6 


1.16 ± 0.02 


7322 


GJ 1111 


08:29:49.35 


+26:46:33.7 


M6.5 


1.59 


7322 


GJ 473b 


12:33:19.1 


+09:01:10 


M7 


1.67 ± 0.01 


7830 


LHS 3003 


14:56:38.31 


-28:09:47.4 


M7 


2.03 


7322 


VR 8 

VD O 


lu.oo.Otj.zy 


D8-9^-zLn 1 

-UO.iO.4U. 1 


M7 


1 .uo 


7^99 


VB 10 


19:16:57.66 


+05:09:00.4 


M8 


2.10 


7322 


Teide 1 


13:05:40.18 


-25:41:06.0 


M8 


1.69 ± 0.03 


9846 


GL 569b 


15:43:03.8 


+26:15:59 


M8.5 


2.29 


7322 


LHS 2924 


14:28:43.33 


+33:10:37.9 


M9 


2.42 


7322 


Roque 33 


03:48:49.0 


+24:20:25 


L0 


1.52 ± 0.06 


9846 


Roque 25 


03:48:30.6 


+22:44:50 


L0 


2.16 ± 0.22 


9846 


Kelu 1 


13:05:40.18 


-25:41:06.0 


L2 


2.79 


7830 


J1228-1547 


12:28:15.23 


-15:47:34.2 


L5 


3.51 


7830 


J0205.5-1159 


02:05:29.40 


-11:59:29.7 


L7 


3.71 ± 0.22 


9846 



a Sp ectra l Types are from IKirkpatrick et al.l (1199 ll ). iHenry et al.l (119941 ). iMartin et al 
j2000h and IKirkpatrick et all Jl999h 
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b Errors for Q are only given for objects with 4 or more spectra 



Table 4. Properties of NGC 1333 objects with spectral types 



Object Name 3 


Spectral Type 


Log(Teff) b 


J-H c 


H 


-K 


K 


A„ d 


Log(L/L ) e 








Q 






Notes 8 


AbKz4 


A /TO _i_ n o 

Me ± U.o 


o.oo 


n ftn 
U.oU 


U 


4/ 


i/.y4 


n a 
U.U 


1 K O 

-l.oz 


n no 
U.Uz 


z 


no 

Uz 


1 


A 

u 


n k 

Uo 


TV /TO O 

Mo.z 


AoKzo 


M4 ± U.O 


o.ol 


n 7n 
U. IV 


A 

u 


fU 


1 £ QQ 

lo.oo 


n v 
U. 1 


-z.ol 


n a o 
U.4z 


1 


9ft 

oO 


1 

± 


n 
U 


nn 

uy 




A QT?00 


Mo ± U.o 


Q /1ft 
O.40 


z.Uo 


1 
1 


10 


1Q OR 

lo.Uo 


lo. f 


i no 
-l.UZ 


n n/i 
-U.U4 


l 


( 


1 

it 


n 
U 


U I 




a qdi nn 


A /TQ i a o 

Mo ± U.o 


o.4U 


1.U1 


A 

u 


1 


1 Q O /I 

lo.Z4 


o.O 


1 no 

-i.yz 


-U.U1 


1 


no 

yo 


1 

± 


n 
U 


n c 
UO 




A CT3 1 1 


M4 ± U.D 


O.OU 


1.10 


A 

u 


c Q 

oo 


14. oy 


4.8 


-z.Uo 


U.UU 


1 


z f 


1 

± 


n 
U 


no 
Uz 




AbKlUo 


Mo ± U.o 


o.4o 


i txft 
1.00 


1 


iy 


1 o ftn 


O K 

o.o 


1 1 

-l.zl 


n on 

u.zy 


1 


ftO 

oz 


1 

± 


A 

U 


H7 

U I 




ASR17 


M6 ± 0.7 


3.44 


1.36 





75 


13.21 


4.8 


-1.37 


0.08 


1 


67 


± 





02 


M7.4 


ASR26 


<M0 ± 0.7 


3.58 


1.59 





85 


15.54 


8.8 


-1.83 


0.18 





90 


± 





01 




ASR8 


M7 ± 1.1 


3.43 


0.63 





33 


12.34 


0.1 


-1.26 


-0.07 


1 


76 


± 





17 




ASR28 


M7 ± 0.8 


3.41 


0.64 





39 


14.18 


0.2 


-2.00 


-0.02 


1 


89 


± 





05 


Binary 


ASR15 


M6 ± 0.7 


3.44 


0.99 





53 


13.49 


3.3 


-1.60 


-0.04 


1 


66 


± 





03 


M7.4 


ASR23 


<M0 ± 0.6 


3.58 


1.48 





76 


13.33 


7.8 


-1.04 


0.16 





89 


± 





07 




MBO20 


M3 ± 0.6 


3.51 


1.75 





88 


10.91 


10.9 


-0.17 


-0.05 


1 


19 


± 





02 




MB078 


M5 ± 0.7 


3.47 


1.78 


1 


18 


13.37 


10.8 


-1.35 


0.16 


1 


45 


± 





05 





a ASR designations are from lAspin et al.l (119941 ); MBO designations are from IWilking et al.l (120041 ) 
b Typical errors are 0.1 

C 2MASS photometry for the objects listed. This was used to derive r^. 
d Typical errors are ~ 1 



e Typical errors are ~ 0.15 
f Typical errors are ~ 0.1 



g Spectral types derived by I Wilking et al.l (120041 ) 
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Table 5. Ratios of intermediate-mass stars to low-mass objects 



Region 


N(0.1-l M )/N(O.O3 - 0.1 M ) 


Pleiades 11 


3.3 +2.2/-1.0 


Chamaeleon b 


5.6 +2.4/-4.5 


Mon R2 C 


5.3 ± 3.35 


Taurus d 


4.0 ± 1.2 


IC 348 e 


3.0 ± 0.5 


Orion f 


3.3 ± 0.6 


NGC 1333 


Ri = 0.14 ± 0.047 


NGC 1333 


R 2 = 0.64 ± 0.38 



a Data from Moraux et all (120031 ) 
b Data from iLuhmanl (]2004al ) 



c Data from lAndersen et al.l (120061 ) 
d Data from ILuhmanl (I2004bl ) 



e Data from iLuhman et al.l (j2003al ) 

f Data from iHillenbrand fc Carpenter! (120001 ) 



